Crossover from dilute-Kondo system to heavy-fermion system 
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Ground state properties of a Kondo lattice model with random configuration of / electrons are in- 
vestigated with a variational Monte Carlo method. We show that the crossover from a dilute-Kondo 
system to a heavy-fermion system occurs when the density of c and / electrons (uc, Uf) become 
comparable, nf < Uc. In the heavy-fermion region, the correlation between / electrons is strong 
and the / electrons themselves greatly contribute to the screening of other /-electron moments. We 
propose that the character of Kondo screening changes from "individual" to "collective" across the 
crossover. 
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Heavy-fermion materials show various interesting phe- 
nomena and have been extensively studied so far. As is 
well known, / electrons with a strong localized character 
play an important role in this system, especially for a 
large effective mass. Origins of this large effective mass 
are considered to be the following factors: hybridization 
between / electrons and other conduction (c) electrons, 
strong Coulomb interaction between / electrons, and a 
certain number of / electrons. In this letter, we focus 
on the third factor, namely, relation between the large 
effective mass and the number of / electrons. 

When some magnetic ions with / electrons are replaced 
by the non-magnetic ions, e.g., Ce^La or Yb— >Lu, the 
system gradually changes from the heavy Fermi liquid 
(dense limit) to the local Fermi liquid (single-impurity 
limit) [E S i. However, the details of this crossover 
are quite different in each compound. For example, 
in Cei_a;Laa:Pb3, the resistivity shows local-Fermi-liquid 
behavior even at x = 0.2 0. It is interesting that the 
single-impurity picture is valid even in a rather "dense" 
region. On the other hand, the situation is rather com- 
plex in Cei_a;Laa; Coins Between the local Fermi liq- 
uid and heavy Fermi liquid, there is an intermediate re- 
gion where a temperature scale characterizing the Fermi 
liquid is greatly suppressed. It is proposed that the spe- 
cific heat, magnetic susceptibility, and resistivity can be 
described by a sum of single-impurity and heavy-fermion 
contributions 0, [Hj- This "two fluid description" indi- 
cates that the system is fluctuating between the local 
Fermi liquid and heavy Fermi liquid in the intermediate 
region. However, the detailed properties of this region 
has not been clarified so far. 

For this problem, several model Hamiltonians, such 
as a periodic Anderson modeWPAM) @, 0, 1 i and a 
Kondo lattice model (KLM) [13, [IH , have been investi- 
gated. A crossover from the dilulte-Kondo system to the 
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heavy-fermion system is well described by these studies. 
However, the character of Kondo screening and the cor- 
relation between / electrons have not been discussed in 
detail. 

In this letter, we study the ground state properties of 
the KLM with random configurations of / electrons, par- 
ticularly focusing on the Kondo screening and the corre- 
lation between / electrons. We show that the crossover 
from the dilute-Kondo system to the heavy-fermion sys- 
tem occurs when the density of c and / electrons {ric, 
nf) become comparable, n/ < ric. In the dilute-Kondo 
region, the / electrons behave independently and the 
single-impurity picture is valid. In the heavy-fermion re- 
gion, the correlation between / electrons becomes strong 
and the / electrons themselves greatly contribute to the 
screening of other /-electron moments. We show the 
change in the character of Kondo screening from "in- 
dividual" to "collective" across the crossover. The vari- 
ational Monte Carlo (VMC) method is used for calcula- 
tion. This method makes it possible to treat the localized 
/ electrons correctly and to take into account the effect of 
randomness beyond the averaged treatment of coherent 
potential approximation. 

We consider the following KLM in a two-dimensional 
square lattice, 

H ^^Ekcl^Cka + J^Si ■ s,. (1) 
k(T ieA 

where Ek is the energy dispersion of c electrons and 
c\.„{cka) is a creation (annihilation) operator of c elec- 
trons with momentum k and spin a. Si and Si 
represent the /- and c-electron spins, namely. Si = 

\ ftf^aa'ha' and S, = \ J^aa' cL^T^^' C<<t' ■ J {> 0) 

denotes the antiferromagnetic exchange coupling between 
them. The / electrons are completely localized and 
treated as localized spins. We call the lattice sites with 
and without / electrons A and B sites, respectively. To 
calculate the physical quantities, we take the average over 
10-40 random configurations of A and B sites. 
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FIG. 1: Schematic phase diagram of the KLM in the strong 
coupling limit as a function of rij. 
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FIG. 2: n/ dependence of a kinetic energy -Bkin and an ex- 
change energy _Bcx for = 0.48 and = 0.84 (J/t = 3.0). 



First, we consider the strong-coupling limit (J/t ^ 1). 
We fix the density of c electrons per site (ric) and change 
the density of / electrons ny (0 < ny: < 1) on the lattice 
with sites, as shown in Fig. [1] In the strong-couphng 
limit, c electrons at an A site forms an on-site singlet 
pair with the localized /-electron spin. When Uf < n^, 
all localized spins form the on-site singlets and remaining 
Ns{nc — nf)c electrons become mobile carriers (called as 
metal 1). On the other hand, when Uf > Uc, c electrons 
are "exhausted" and Ns{nf — ric) unscreened local- 
ized spins exist. In this case, when a c electron moves to 
a site with an unscreened localized spin (unscreened A 
site), the original on-site singlet disappears and another 
one is newly created. Thus the hopping of c electrons 
is described as a motion of on-site singlet among the A 
sites. It will cover all the A sites and will dynamically 
screen all the localized spins. We can say that this screen- 
ing is "collective" in the sense that / electrons effectively 
participate the screening of other /-electron moments. 
The motion of on-site singlets can be regarded as a mo- 
tion of mobile holes at the unscreened A sites and it is 
a different picture with metal 1 (metal 2). In the dense 
limit [uf = 1), this system is identical to the U = oo 
Hubbard model with a reduced (half) band width p^ . 
Thus the metal 2 is considered to be a simple model for 
heavy-Fermi liquid. The metal 1 and metal 2 have dif- 
ferent characters and separated by the special density of 
Uf — fic (see Fig. [TJ , where the system has a gap of or- 
der J/t and becomes insulator. For a finite value of J/t, 
this simple insulating state is not expected but it is pos- 
sible that the character of metallic state changes around 

For finite J/t, we use the VMC method to calculate 
the physical quantities of the ground state. We prepare 
the Gutzwiller-type trial wave function jvl/) for a trial 
state [3]. It consists of a projection operator Pj and a 
one-body part 1$) that is obtained from the solution of 
a certain one-body Hamiltonian. The specific form is as 



follows, 

I*) = Pf 1$) , Pf^U [4(1 - Ml) + 4(1 - 4)" • 

(2) 

Pf keeps the /-electron number of each A site exactly 
one. Since the / electrons must be treated as localized 
spins in the KLM, this constraint is necessary for the 
accurate calculation The one-body part |$) is ob- 
tained by diagonalizing the U = PAM in a real-space 
representation. The explicit form is as follows, 

{i,j)cr i,o" i,o" 

(3) 

Here, Uj^ are transfer integrals from site j to i with 
spin (T. Via and Ef. are "effective" hybridization and 
/-electron energy level for site i with spin a. To avoid 
the occupancy of / electrons at B sites, we set Ef^^ = oo 
(or a sufficiently large value in the actual calculation) for 
I S B. tija, Via and Ef_^^ {i £ A) are variational pa- 
rameters and optimized so as to minimize the variational 
energy. Here, we set Uja = i, Via = V and Ef.^ = Ef 
for simplicity. Then we calculate physical quantities with 
the optimized wave function. To obtain the expectation 
values, we take 10-40 patterns of random configurations 
and average them. Although a large number of different 
configurations are necessary, this calculation procedure 
will take account of the effect of randomness more pre- 
cisely than the coherent potential approximation that is 
often used for disordered systems. 

We perform a VMC calculation for a square lattice 
with Ek = —2t{coskx + cos ky). i = 1 is set to be a unit 
of energy in the following. The size dependence of results 
have been checked from 10x10 to 16x16. In the system 
with 10x10 sites, we have found some finite-size effects 
at n/ < 0.16 [l^ but they disappear in the system with 
more than 14x14 sites. Since the 10x10 system is not 
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FIG. 4: n/ dependence of C^'" and C^~^ for (a) = 0.48 
and (b) — 0.84 (J/t = 3.0). Dotted lines represent the 
strong-coupling limit. Arrows indicate the crossover points. 



FIG. S: Uf dependence of strength of on-site screening S for 
Uc = 0.24, 0.48 and 0.84 {J/t = 3.0). Dotted lines repre- 
sent Uf 
points. 



0.24, 0.48 and 0.84. Arrows indicate the crossover 



sufficient for describing the dilute region, we show the 
results of 14x14 in the following. 

First, Fig. [5] shows the Uf dependence of kinetic en- 
ergy fi'kin and exchange energy £^ex that are expectation 
values of the first and second term of right-hand side of 
Eq. ([1]). Uf ^ corresponds to the dilute limit and 
n/ = 1 corresponds to the usual KLM. For the case with 
Tic — 0.48 (Fig. [2)Ja)), i?ox shows almost linear behavior 
until Hf ~ 0.40, namely, the exchange energy per / elec- 
tron (Ecx/Nf) is almost unchanged even with increasing 
the number of / electrons. It means that the additional / 
electrons do not affect the screening of other / electrons 
so much. We can say that the / electrons in this region 
behave independently and the single-impurity picture is 
valid until n/ ~ 0.40. Because of the finite value of J/t, 
this boundary is slightly different from that expected in 
the strong-coupling limit, 0.48. The effect of finite J/t is 
discussed later. For the case with ric = 0.84 (Fig. [2I^b)), 
as can be seen by the abrupt changes in -Ekin and -Eex, 
the single-impurity picture is valid until ?i/ — 0.84. It is 
almost the same with the strong-coupling limit. 

Next, we calculate the following quantity. 



S 



JNf 




(4) 



rnd 



where (• • • )^^^^ denotes the random average. This quan- 
tity is always positive and corresponds to the strength 
of on-site screening per / electron. We fix the value of 
ric as 0.24, 0.48 and 0.84 and show the Uf dependence 



of S for J/t = 3.0 in Fig. H For = 0.84, as n/ in- 
creases, S gradually increases and rapidly decreases at 
Uf ~ He = 0.84. This crossover at n/ ~ nc is con- 
sistent with the strong-coupling limit, and thus corre- 
sponds to the crossover from the dilute-Kondo system to 
the heavy-Fermion system. As shown later, the character 
of screening changes from individual to collective at this 
density. For ric = 0.48, on the other hand, the crossover 
occurs at n/ ~ 0.40, slightly away from Uc- We consider 
that this is because the screening clouds have a finite 
spread for a finite value of J/t. In the strong-coupling 
limit, the screening clouds are completely localized at 
each site and have no overlap. In this case, the bound- 
ary between the dilute-Kondo system (individual screen- 
ing) and heavy-fcrmion system (collective screening) is 
exactly at ~ Uc- However, the screening can be collec- 
tive if the neighboring screening clouds overlap and inter- 
act with each other, even though nj < nc- For Uc = 0.48, 
the spread of screening cloud will exceed the lattice spac- 
ing and as a result, the region of heavy- fermion system 
extends to n/ < Uc- For ric = 0.24, the change of S be- 
comes more gradual and we cannot determine the clear 
crossover point from these data. We conclude that the 
more Uc decreases, the more the screening cloud spreads 
and the crossover becomes ambiguous. 

To observe the character of screening more clearly, we 
study the correlation between spins. For each /-electron 
spin Si, the following relation is satisfied. 



f-c 



C: 



■f-f 



+ - 



(5) 



because Si forms singlet pairing with the rest of all spins 
in a paramagnetic state. We calculate C-^"^ = (c/~'^) 
and C^^^ — {C(~^), where (• • • ) denotes an average 
over both sites and random configurations. C^~'^ {C^~^) 
indicates how much the c (/) electrons are concerning 
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with the screening. Figure [4] shows the n f dependence 
of these quantities for ric — 0.48 and 0.84. In the single- 
impurity hmit {nf 0), screening is performed only 
by c electrons and C^'" 3/4 (0.75) and C^'f 0. 
As n/ increases, / electrons begins to participate the 
screening and C^^^ increases while C-^~'^ decreases. At 
a certain point around Uf ^ Uc (indicated by arrows in 
Fig. HI), an abrupt change is observed. It corresponds to 
the crossover point shown in Fig. [3] and represents the 
crossover from the dilute-Kondo system to the heavy- 
fcrmion system. In the heavy-fermion system, / elec- 
trons themselves are greatly participating the screening 
of other /-electron moments and therefore the correlation 
between / electrons is much stronger than in the dilute- 
Kondo system. We consider that this strong correlation 
between / electrons is an origin of large effective mass in 
the heavy-fermion system. 

Here, we compare our result with other theoretical 
works. Mean-field calculation by Kaul and Vojta for 
20x20 square lattice has shown the n/ dependence 
of T* = (6^) /D, where 2br = J {c\.fr) denotes auxihary 
fields and Z) denotes a half band width. T* is almost 
the same quantity with S in our result and also shows 
the crossover around Uf ~ Uc- Burdin and Fulde have 
calculated Tq/Tk using a coherent potential approxima- 
tion and a dynamical mean- field theory approach TT], 
where Tq denotes a coherence (Fermi liquid) tempera- 
ture scale and Tk denotes a Kondo temperature. Tq/Tk 
is an important ratio and related to the effective mass 
as Tq/Tk = "^sikm/'^klM' where SIKM denotes single- 
impurity Kondo model (n/ — > limit of KLM). Their 
result have shown the rapid decrease of Tq/Tk, namely, 
the enhancement of effective mass around — Uc- It 
is consistent with our result showing that the crossover 
from the dilute-Kondo system to the heavy-fermion sys- 
tem occurs around Uf < ric- The enhancement of ef- 
fective mass for > Uc is also shown in the PAM by 
Pruschke et al 20]. 

Finally, we refer to the relation between our result and 
the experiments. Ybi_2;Lu2;Al3 shows clear deviation 
from the single-impurity model up to a; 0.7 («/ ~ 0.3) 
and the heavy-fermion picture is valid in this region [13] . 
The estimation of ric ~ 0.5 in this compound [l^ does 
not contradict our result. The systematic experiment of 
YbXCu4 has shown that the compounds with large ric 



(X=Ag, Tl) are well described by the single-impurity 
model, while the compounds with small ric (^=Mg, Zn, 
Cs) arc not [19]. This result is reasonable at least quali- 
tatively if the former correspond to the dilute-Kondo sys- 
tem {rif < He) and the latter correspond to the heavy- 
fermion system (rif > ric). In Cei„a;La2:CoIn5, an in- 
termediate region with both the single-impurity and the 
heavy-fermion contributions are observed ,4]. This indi- 
cates that the crossover from the dilute-Kondo system to 
the heavy-fermion system is rather ambiguous and the 
system is fluctuating between them. In our result, the 
ambiguous crossover appears in the case of small J/t (not 
shown in this paper) and small ric because the screening 
cloud is fairly extended in this case. We expect a smaller 
J/t and/or a smaller ric for Cei-^Lax Coins, compared to 
the materials with clear crossover, such as Cei„a;LaxPb3. 

In this paper, we have studied the ground state prop- 
erties of the KLM with random configuration of / elec- 
trons. The VMC method in real space is used for cal- 
culations. We have found that the crossover from the 
dilute-Kondo system to the heavy-fermion system occurs 
around rif <nc. In the dilute-Kondo region, the / elec- 
trons behave independently and the single-impurity pic- 
ture is valid. In the heavy-fermion region, on the other 
hand, the correlation between / electrons is strong and 
the / electrons themselves participate the screening of 
other /-electron moments. It is considered to be an ori- 
gin of large effective mass. We propose that the screen- 
ing character changes from "individual" to "collective" 
across the crossover. 
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